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It has been very recently reported (McKeen, Pospelov, arXive:1011.3046) that the parameter 
region suggested for an explanation of the neutrino oscillation results from the LSND, KARMEN 
and MiniBooNE experiments in terms of the production and radiative decay of a heavy neutrino (vh) 
can be ruled out based on the measurements of the radiative muon capture (RMC) on hydrogen. We 
calculate limits on mixing strength between the Vh and by using results of this experiment, and 
find out that they essentially disagree with the reported bounds. For the vu with mass 60 MeV our 
limit is worse by a factor of 5, while for 90 MeV it is worse by about two orders of magnitude. We 
also noticed the wrong behavior of the reported limit curve in the low mass region. The importance 
of accurate Monte Carlo simulations of the vu signal in the RMC experiment is stressed. Our 
conclusion is that the whole LSND-MiniBooNE parameter region cannot be ruled out by the RMC 
measurements. 



PACS numbers: 

In the recent work [![, see also @, it has been shown 
that the neutrino oscillation results from the LSND, 
KARMEN and MiniBooNE experiments could be ex- 
plained by the existence of a heavy sterile neutrino (vh), 
assuming that it is created by mixing in v^ neutral- 
current interactions and decays radiatively into a photon 
and a light neutrino. The Vh's could be Dirac or Ma- 
jorana type and decay dominantly into Vh —> jv in the 
detector target if, for example, there is a non-zero transi- 
tion magnetic moment between the Vfi an d the v. Com- 
bined analysis of the energy and angular distributions of 
the LSND and MiniBooNe excess events suggests that 
the Vh mass m Vh , the mixing strength 1/7^ | 2 and the 
lifetime r VK , are in the ranges: 

40 < m Uh < 80 MeV, 10~ 3 < \U^ h \ 2 < 10~ 2 , 

10~ n < r Vh < KT 9 a, (1) 

respectively. 

The mixing l?/^ 2 for the mass range of ([T]) would 
result in the Vh emission in the ordinary muon capture 
(OMC) on nuclei pT A -> v^A' @. The OMC rate of the 
heavy neutrino production can be estimated as 

r Uh =ToMcK h \ 2 p(m Uh )/p(0) (2) 

where Yomc is the OMC rate and p{m Vh ), p(0) are the 
phase space factors for emission of heavy neutrino and 
massless, unmixed neutrinos, respectively. 
Recently, it has been noticed by McKeen and Pospelov 
P[ that the parameter space of ((J) can be probed by us- 
ing the results of the experiment on the radiative muon 
capture (RMC) rate on hydrogen Q. In this experiment 
muons were stopped in a liquid hydrogen target. Pho- 
tons from the reaction pp — > v^jn were converted in a 
Pb layer surrounding the target into e + e~ pairs , whose 
momenta were measured by a magnetic spectrometer. If 
the Vh exists, it would be produced through mixing in 
the OMC process and after the prompt Vh — > "fv decay 



in the target and subsequent decay photon conversion re- 
sult in a final state identical to the one from the RMC 
reaction. The number of RMC photons observed in the 
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FIG. 1: Energy distributions of photons from the vu — ¥ 
decay of the Dirac heavy heutrino with a = — 1 [l[ with energy 
of 100 MeV created through the mixing in the muon capture 
H v p — > Vhti on proton and calculated for Vh masses of 60 and 
80 MeV. The arrow indicate the threshold of 60 MeV used in 
the analysis of RMC data 0|. 

experiment is calculated by using Eq.(38) from @: 

n 7 = n stop R 7 KA a (3) 

where n s t op {— 3.4 ■ 10 12 ) is the total number of muons 
stopped in the target, i? 7 = (2.1±0.21)-1(T 8 is the RMC 
branching fraction, K — 0.6 is an effciency factor, and A a 
is the absolute photon spectrometer acceptance which 
include geometrical acceptance, conversion probability, 
e+e~ reconstruction efficiency etc.. Using n 7 = 279 ±26 
events found in the experiment, gives averaged over the 



2 



RMS photon spectrum value A a = 6.5 ■ 10 -3 . We will use 
these numbers for the cross-check, but immediately note 
that the A a value is about two times bigger than that 
estimated for the Uh — > "(v spectrum. To be conservative 
we will keep this "safety" factor just in mind. 
The energy distribution of photons from the Vh — > jv de- 
cay calculated for different masses of the Vh with the total 
energy of 100 MeV, which is expected from the OMC on 
hydrogen, is shown in FigJT]The upper limit on the excess 
number of photons from the v^, ~ > 71/ decay in the RMC 
experiment can be estimated at 2tr level for E~ > 60 
MeV as: 

An 7 = n 7 — 7j < 140 events (4) 

where no ~ 180 is the number of events that are expected 
from the theory for the ratio of the pseudoscalar and 
axial- vector form factors g p /g a = 0.69 and ortho to para 
transition rate of the (/ip) atom equal A = 4.1 x 10 s _1 

1- 

The excess n s of events from the Vh — s> jis decays can be 
calculated as 

n s = n stop R mc\U phi 2 f P h s B(v h ^v)P Vh ^ lu f 1 KA a 

(5) 

where n stop , K, A a are as in (O, i? Mc(L5 ■ 10 3 ) 
is the OMC branchning fraction, B(vh — > 7^)(~ 1) 
is the branching fraction of the Vh — > 7^ decay, 
fphs = p{ m v h ) I 'p(O) is the ratio of phase space factors, 
^ h ^ 7 i/(0.15) is the probability to decay in the target for 
heavy neutrino lifetime r = 10~ 9 s as in [4], / 7 is the frac- 
tion of photons with energy E 1 > 60 MeV (0.074 and 0.05 
for 60 and 80 MeV Vh , respectively) . Taking into account 
above values results in the 2a limit l^fil 2 < 1.4 x 10~ 3 
for the Vh mass of 60 MeV, which is worse than the limit 
of [H by a factor of 5, see Fig. 3 in For the mass of 80 
MeV the limit is < 4 x 10~ 3 and the disagreement 

is about factor 10, but for 90 MeV it is two orders of mag- 
nitude, \U^ h \ 2 < 5.1 x 10~ 2 vs |[/ M | 2 < 5 x 10- 4 reported 
in For the mass range m Vh < 20 MeV the behavior 
of the limit curve shown in Fig. 3 Qj and its extrapola- 
tion to zero mass is incorrect , because the probability 
Pv h ^t~fv — > with m Vh — > 0. 

Let us note that work of Ref.[l[ has been attempted to 
explain an excess of events observed by the LSND in the 
energy range 20 < E 1 < 60 MeV. The idea is to in- 
troduce a new particle with such decay properties that 
would allow to keep the number of events with the energy 
deposition above 60 MeV close to zero, in agreement with 
the LSND observations. This is achieved by suggesting 
that the Vh decays radiatively as a Dirac particle with 



the photon asymmetry parameter equal to a = — 1. The 
emission of photons preferably backward with respect to 
the Vh direction of move makes their energy spectrum 
much softer compare e.g. to the one form the isotropic 
distribution (probably used in the analysis of Ref.[4|). In 
contrary the limits of Ref . [I[ are extracted from the data 
of the RMC experiment for the energy region E 1 > 60 
MeV and even for the lower Vh energy than in the LSND 
case. The fraction of expected Vh — > jv events in this 
energy range is small, and hence the limits obtained are 
quite sensitive to the details of the experimental analy- 
sis, in particular for the Vh masses above 60 MeV. Let us 
give an example of such sensitivity related to the photon 
response function in the RMC experiment. It is deter- 
mined by generating Monte Carlo photons sampled from 
the known 7r~p reaction spectrum, whose origins in the 
target are identical with the pion stopping distribution. 
The shape of this function is important to extract cor- 
rect number of RMC events. The function is essentially 
asymmetric ( see Fig. 10), and shifts the whole energy 
spectrum to the low energy region. The peak value of 
129 MeV is shifted by 5-7 MeV but the low energy tail 
has energies up to 50% of the initial photon energy. If wc 
shift spectra shown in Fig.l just by, say 4 MeV, the de- 
crease in the number of observed Vh — > jv events would 
result in limits |[/ M ,,| 2 < 2.1 • 10~ 3 and |[/ M | 2 < 0.01 
for the 60 and 80 MeV Uh, respectively, instead of given 
above. Therefore, accurate Monte Carlo simulations of 
the Vh production and decay sequence events in the ex- 
periment [5J and propagation of theses events through 
the spectrometer taking into account the response func- 
tion is important to extract reliable limits on l/J^I 2 - 
The Primakoff mechanism of the Vh production vA — > 
VhA followed by the decay Vh — > "fv through the tran- 
sition magnetic moment \itr between the Vh and light v 
with the subsequent Vh — >• jv decay has been originally 
considered in [6j and has been used to extract limits on 
Ht r from the NOMAD data 0] ■ In order to use this mech- 
anism for interpretation of the LSND v^ data [i| , one has 
to make an additional assumption that the light v is a 
component of the v^. Although this is possible, it re- 
quires introduction of additional unknown parameters. 

In summary, we do not find arguments of McKeen 
and Pospelov [J] convincing enough to exclude the whole 
LSND-MiniBooNE parameter region, at least for Dirac 
scenario of [l(. In contrary, their work enhances motiva- 
tion for a more sensitive direct search for the Vh — > 7^ 
decay in a dedicated experiment. 

I would like to thank D.S. Gorbunov for discussion and 
pointing out Ref. Q to me. 
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